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ABSTRACT 

From the results of numerical scattering experiments and simulations of a massive 
black hole binary in spherically symmetric and shallow cores it has been deduced that 
most likely the shrinking process stalls due to loss-cone depletion before emission of 
gravitational radiation becomes important. Here we follow a different approach and 
focus on the population of stars which is bound to the binary and so far has not 
received much attention. With simple assumptions which should not be sensitive to 
initial conditions we derive a lower limit for the mass of stars which needs to be ejected 
by the binary in order to coalesce. We also compute this mass in dependence on the 
steepness of the density profile according to which the stars are distributed. Our results 
are not as pessimistic as earlier conclusions and actually suggest that the BHs merge. 

Key words: black hole physics - galaxies: evolution - galaxies: interactions - galaxies: 
kinematics and dynamics 
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1 INTRODUCTION 

The existence of a massive black hole binary (BHB) is a nat- 
ural consequence of the assumptions that galaxies harbour 
a massive black hole (BH) in their center and that galax- 
ies merge with each other. The evolution of the binary can 
be split up in roughly three phases: First both cores with 
the BHs in their center spiral inwards to each other due 
to dynamical friction. When the BHs bind to each other 
and form a hard binary on the parsec scale, the semimajor 
axis continues to decay due to slingshot ejection of stars. 
Eventually, in the third phase, emission of gravitational ra- 
diation dominates further shrinking until the BHs coalesce 
(Begelman et al. 1980) . It is still a matter of debate whether 
the BHs coalesce or the shrinking stalls before entering the 
final phase because no stars are left to interact with the 
binary (i.e. loss-cone depletion). However, the origin of X- 
and Z-shaped radio galaxies is probably best explained by 
the coalescence of the two BHs (Rottmann 2001; Zier & 
Biermann 2002; Gopal-Krishna, et al. 2003; Zier 2005) and 
the observed number of these sources is in agreement with 
the merging rate derived for radio galaxies (Merritt & Ekers 
2002). Numerical scattering experiments showed that the 
BHs can merge on scales of ~ 10 8 ~ 9 yr if the loss-cone is 
always full (Quinlan 1996; Quinlan & Hernquist 1997; Zier 
k, Biermann 2001; Milosavljevic & Merritt 2001). However, 
for a flat core most authors estimate that the loss-cone gets 
depleted long before the binary can enter the third phase. 
Further hardening then depends on the rate at which the 
loss-cone is refilled, most probably by two-body scattering 
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of stars. This implicates a time scale for the separation of the 
BHs to shrink to a distance where gravitational radiation be- 
comes dominant which might exceed a Hubble time so that 
the binary basically stalls. The conclusions are essentially 
the same: though the results show that the BHs coalesce 
it is argued that in real galaxies the loss-cone is depleted 
before gravitational radiation dominates further shrinking. 
These conclusions are based on the assumption that at the 
time the binary becomes hard the density profile of the cusp 
is flat (i.e. with a power-law index of at most 2, usually 1 or 
even 0). This seems to have been confirmed by recent simu- 
lations of a binary in a constant density core (Berczik et al. 
2005). As the authors say themselves the choice for the mass 
of the binary is unrealistically large compared to that of the 
galaxy and the initial conditions are quite unlikely, because 
the BHs were introduced symmetrically about the center of 
the galaxy instead of arriving there due to the evolution of 
the merger. 

The initial conditions in these numerical experiments 
are idealised and cover only a small fraction of the param- 
eter space. They are based on spherically symmetric pro- 
files which have been derived from observations of elliptical 
galaxies. These profiles correspond to the central cluster be- 
fore or after the merger when mass has already been ejected 
from the loss-cone and shifted from smaller to larger radii, 
resulting in flatter profiles. However, they do not match the 
central profiles during the merger. Also ignored are individ- 
ual spins of both galaxies, which might stabilize the cluster 
against tidal disruption by the other BH beyond a hard bi- 
nary, depending on the magnitude of the spins and their 
orientation relative to each other and relative to the orbital 
angular momentum of the merging galaxies. These param- 
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eters have a strong influence on the merger itself and the 
morphology of the remnant, as has been shown by Toomrc 
& Toomre (1972). When galaxies collide energy will be dissi- 
pated and angular momentum redistributed with some frac- 
tions compensating each other. Large amounts of mass will 
move on highly eccentric orbits (Rauch & Tremaine 1996) 
in a potential that is strongly non-spherically symmetric so 
that the angular momentum of a single particle is not con- 
served and matter with low angular momentum piles up in 
the center. The density of both cores will be increased con- 
siderably before they merge and the binary becomes hard 
(Barnes & Hernquist 1996). Each of the BHs will carry a 
stellar cusp as massive as the BH and therefore a stellar 
mass comparable to that of the binary will be concentrated 
in the central core once the cusps merge. Hence we expect 
that at the time the BHs form a hard binary, the surround- 
ing density distribution will be much more compact with a 
steeper profile than in non-interacting galaxies, i.e. unlike 
the initial conditions used so far in numerical experiments. 
It is the fraction of this cusp which is bound to the binary 
to which our analyses applies. 

Here we will study in a simple approach how compact 
such a cusp has to be and wich profile is required so that 
the binary can shrink to the third phase and the BHs even- 
tually coalesce. Whereas numerical experiments were more 
concerned with unbound stars scattered off the binary, we 
will focus on stars which are bound by the BHs. Therefore 
our results will be less sensitive to initial conditions. 



2 PRELIMINARIES 

Being interested in the 2nd phase of the merger we assume 
that the BHs have bound to each other and are moving on 
Keplerian orbits. The origin is the center of mass of the bi- 
nary and we define the mass ratio q = m2/m\ < 1. The total 
and reduced mass are M12 = mi + m2 and fi — mim.2/Mi2, 
respectively. With a we denote the semimajor axis of the bi- 
nary, i.e. of the orbit of the reduced mass. For circular orbits 
this is equal to the separation of the BHs. Thus the energy 
of the binary can be written as 
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and the relative velocity between the BHs is 



GM 12 



(1) 



(2) 



which corresponds to the velocity of the reduced mass if it 
moves on circular orbits. The definition for the semimajor 
axis ah where the binary becomes hard is not unique and 
changes in the literature (e.g. Heggie 1975; Hills 1975; Hut 
1983; Quinlan 1996; Milosavljevic & Merritt 2001). Some 
definitions are derived from the results of numerical experi- 
ments and are more phenomenological. Usually they are all 
based on comparing the velocity dispersion of the cluster 
(which itself is not so easy to define in an ongoing merger 
when the core is far from being relaxed) with the velocity of 
a component of the binary. For major mergers, i.e. large q, 
they all yield similar results for ah, which is found to be on 
the parsec scale for Mi ~ 10 8 Mq. 

We think the transition from the first to the second 
phase is best defined by the distance between the BHs when 



the stars are moving in the potential of both BHs, not just 
one, so that the binary and the star can be treated as a re- 
stricted three body problem. This transition from dynamical 
friction to slingshot ejection as dominating processes for the 
decay of the binary is smooth. Because we do not know an 
exact definition for a h we instead scale it to the semimajor 
axis a g at the end of the second phase, which is more clearly 
defined. An upper limit for this transition is if it still takes 
a Hubble time for the BHs to merge completely from a = a g 
due to emission of gravitational radiation. For circular orbits 
this is (Peters 1964) 
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In previous publications the ratio of the semimajor axes 
where the transitions between the phases occur, i.e. r\ = 
ah/a g , has been found to be in the range of 20-100. Ap- 
plied to Eq. (3) this yields ah to be in the range 1-7 pc, in 
agreement with a h being on the parsec scale. Because this 
ratio does not seem to depend sensitively on the initial con- 
ditions and is quite robust we will use it in the following to 
compute the location of a^. Note that ah and a g scale differ- 
ently with the mass of the BHs and the velocity dispersion 
of the cluster. Consequently the ratio 77 will be a function of 
these quantities. Because we are considering a range of fixed 
values for rj, this does not affect our following analysis and 
might become important only if systems with very different 
BH masses and velocity dispersions are considered. 

For bound stars it makes only a small difference whether 
the stellar distribution and potential is spherical, axisym- 
metric or triaxial. For simplicity we consider a spherical den- 
sity distribution which in the inner region follows a power- 
law, p = po(r/ro)~ 7 . With the total mass of this cluster 
(M c ) being distributed between the inner and the cluster 
radius, n and r c respectively, the mass distribution is 



„ 3 "7 



7/3 



(4) 



ln(r c /n) ' ' 

To avoid a singularity at r — we cut off the distribution at 
the inner radius r\. The circular velocity is v^ ilc — r d&/dr — 
GM(r)/r. A star with this velocity is bound to a BH with 
mass M12 at the center of the cluster if the velocity is less 
than the escape velocity v csc = y/2GM\2/r. Using the cir- 
cular velocity as the typical velocity (i.e. the velocity at the 
maximum of an isothermal sphere with a Maxwellian veloc- 
ity distribution, same as v^-times the velocity dispersion of 
this distribution) we obtain from « c i rc < v CBC the relation 
M(r) < 2M12. Thus a mass of about 2M12 of the cluster is 
bound to the binary. We expect a large fraction of this mass 
to be in the loss-cone. 



3 REQUIRED MASS AND ITS DISTRIBUTION 

In this section we will derive limits for the mass which is re- 
quired to be ejected and how steep the mass distribution has 
to be so that the BHs can merge. The inner region where the 
stars are bound to the binary is dominated by the potential 
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of the BHs. Approximating the potential of the binary to 
first order with a point potential of the mass M12 located at 
the cluster's center introduces only minor deviations whith 
a maximum of a factor less than 2 for q — 1. Before a star 
becomes ejected via the slingshot mechanism its binding en- 
ergy is about 

, GMurn, 
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where r_ denotes the pericenter and e < 1 the eccentricity 
of the orbit. In this expression we neglected the potential of 
the cluster itself, whose mass amounts up to twice the mass 
of the binary, as we have shown in the previous section. The 
influence of the potential of the cluster could be explored 
in self-consistent calculations using a potential generated by 
the cluster until its mass drops below that of the binary. 
We do not expect the cluster potential to change the basic 
results obtained in the present work and leave these calcu- 
lations to a subsequent paper. The final energy of the star 
after its ejection, E, t { , is zero or more. Scaling it to the initial 
energy with e = we can write it as E*j = K GMi2m*/2r_, 
with k being the scaling factor. Independent of the density 
profile Quinlan (1996) finds that the dominant contribution 
to the hardening of the binary comes from stars whose clos- 
est approach to both BHs is about the semimajor axis a. 
Hence we replace r_ with a so that the energy change of the 
star E t j — E f> i can be written as 
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with k being defined by the last equality. From scattering ex- 
periments it follows that k w (3/2) 2 m2/Mi2 (Quinlan 1996) 
or ~ fi/Mi2 (Saslaw et al. 1974). Thus for circular orbits 
(e = 0) we find k in the range 1 < k < 2, or as large as 
3.2 according to Yu (2002). Even for highly eccentric orbits 
(e ~ 0.7) the energy change will be of the order of one. 

3.1 The ejected mass 

When a star is ejected it extracts the amount of energy 
given in Eq. (6) from the binary. In the limit m* <ti ra 2 
we can replace m* with dm and using Eq. (1) we can write 
dEbiu — AE f as 



This is only slightly more than M12 for q = k = 1 and 
77 = 100 and therefore less than the mass bound to the 
binary. In deriving this expression we assumed that the stars 
are moving on orbits with pericenters which are equal to the 
current semimajor axis of the binary. However, stars interior 
to the binary will at least be disturbed, if not ejected, once 
the binary has become hard. These stars are more tightly 
bound in the potential of the BHs and have to gain energy 



on the expense of the binary in order to be shifted to orbits 
with a as pericenter. This issue will be explored in detail 
in the next section which shows that the mass obtained in 
Eq. (8) is sufficient for a merger if it is distributed with a 
power-law index 7 = 3. 

To define the mass ejection rate Quinlan (1996) intro- 
duced a similar expression as Eq. (7) and used as mass scal- 
ing factor M12 without further justification. This is plausible 
for large q because in this case it can be expected that the 
binary has to eject a mass of about its own. This was con- 
firmed by his results which are essentially in agreement with 
ours. However, in the limit m,2 — > m» it does not seem to 
be plausible that a total mass of about M12 is required. The 
energy of the binary is that of the reduced mass orbiting in 
the fixed potential of the point mass M12 with an angular 
momentum of L bin — fiav^. Therefore a mass of about fi, 
as derived in the above equation, seems to be more natural. 
In order to explain mass deficits in galaxy cores Milosavl- 
jevic et al. (2002) find that a scaling factor of 7712 instead of 
M12 matches the ejected mass better. This dependency on 
q is very similar to that in our expression, supporting our 
approach of focussing on the bound population of stars. Be- 
cause we did not make use of any further assumptions in the 
derivation of Eq. (8) it should best describe the true ejected 



3.2 The distribution of the ejected mass 

According to Milosavljevic & Merritt (2001) the results of 
their simulations indicate that soon after the binary becomes 
hard the loss-cone becomes depleted in little more than the 
local crossing time of the cusp. If we assume that all the 
mass of Eq. (8) is ejected instantaneously at a h , i.e. that the 
density is distributed according to a Dirac-delta distribu- 
tion that is non- vanishing at r — ah, the energy change due 
to ejection of this mass is Ai5 e j = (GMi2/u/2oh) ln(oh/a g ) 
(Eq. [6]). Equating this with the change in energy of the 
binary in Eq. (1) if the semimajor axis shrinks from ah to 
the final semimajor axis af 
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giving an expression for the shrinking of the binary da due a{ 
to the ejection of a mass dm. Stars are ejected all the time 
as long as the loss-cone is not depleted and so the binary 
hardens continuously. Integration from ci g to ah yields the 
mass that has to be ejected by the binary in order to enter 
the last phase and we obtain 

nici= k ln ~- (8) 



and solving for their ratio we obtain 
h = 1 + ln(a h /a g ), 



(9) 



(10) 



which is only about 5.6, 4.9 and 4.0 for ah/a g = 100, 50 
and 20 respectively. However, stars at radii smaller than 
ah are also bound deeper in the potential of the binary. 
Before the above derivation for at applies the BHs have to 
transfer sufficient energy to these stars so that they move on 
orbits with a radius of ah. In the point potential the binding 
energy of a mass element dm — 4nr 2 p(r) dr in a distance r 
is dE = GM\2dm/2r. The integration over the cusp's profile 
between a g and ah gives the binding energy of the cusp 
in this potential and corresponds to the energy which the 
binary loses when it ejects these stars. This negelects the 
potential of the cusp itself, and because the final energies 
of the stars are probably more than zero this is actually 
a lower limit for the energy loss of the binary. Using the 
density profile p — po(?Vro)~ 7 , and with the help of Eq. (4), 
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Figure 1. The final semimajor axis in units of a g as function of 
the exponent 7. Solid, dashed and dotted lines show the distribu- 
tion for different ratios rj = a^/a g . For 7 < 4 they split up in a 
lower (A = 1) and upper (A = 10) branch. 



Figure 2. The mass distributed between a g and which is 
required to allow the BHs to shrink to a g in units of m c j as 
function of the exponent 7. (k = 1). 



the integration yields 



AE = 



GM 2 iM c 



—^7- ln(a h /a g ) 



a h 
2-7 



3-7 _ 

2^7 



r 

2- 7 



7 = 2 
7 = 3 

else. 



(11) 



We assume that the lower limit for the mass which is re- 
quired to allow the BHs to merge, found in Eq. (8), is the 
cluster mass M c which is distributed between n and r c ac- 
cording to the above power law. If the binary ejects the 
stars from the region a g < r < a h we can compute the final 
semimajor axis of the binary in the same way as above, i.e. 
equating Eqs. (9) and (11), and find: 
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Here we have defined the following ratios: rj — <2h/a g > 1 
as before, £ = r c /n > 1 and A = r c /ah > 1. If we allow 
n < a g , so that the cluster extends to radii smaller than 
a g , the fraction of mass between these two radii is not in- 
cluded in the ejected mass since we used a g as lower limit in 
the integration (Eq. [11]). The steeper the distribution, the 
larger is this mass fraction, and hence the less the binary will 
shrink so that a,f/a s increases with 7, if 7 exceeds a certain 
value. However, these solutions are unphysical and will not 
be considered. Therefore, setting n = a g , we can substitute 
( — \tj. In Fig. 1 we plotted the result for k — 1 as ratio 
<Zf/a g = r/ai/ah in dependence on the exponent 7. As ex- 
pected, the steeper the cusp the closer the binary shrinks to 
a g . In the range 7 < 4 the solutions split into a lower (A = 1) 



and upper (A = 10) branch. In case of A = 1 the mass of 
Eq. (8) is distributed between a g and a h , which coincide 
with n and r c respectively. For a uniform density distribu- 
tion (7 = 0) the final semimajor axis is smaller by a factor 
of about 1.4 than for the case of the Dirac-delta distribution 
which yields a f /a g w 17.8, 10.2 and 5.0 for rj = 100, 50 and 
20 respectively, see Eq. (10). The decrease of the final semi- 
major axis with increasing 7 is steeper for large r\. If we write 
in Eq. (7) dm — 4irr p(r)dr and solve for the density, we ob- 
tain p(r) = fi/4-Kkr 3 , i.e. the power-law for which m c j is just 
enough to allow the BHs to merge (of/o g = 1 at 7 = 3). For 
a distribution as steep or steeper than p oc r -3 the ejected 
mass is bound deeply enough in the potential of the binary 
prior to its ejection for all ratios rj in order to allow the 
BHs to shrink to the distance where gravitational radiation 
begins to dominate the further decay. Hence for sufficiently 
steep and compact distributions the BHs coalesce. 

If A = 10 the mass m c j is distributed in a larger sphere 
with an outer radius 10 times as large as a h . The larger A 
is the closer the ratio Of/a g approaches r\ at 7 = because 
less mass is available for ejection at radii < ah, and the less 
the binary will shrink (upper branches in Fig. 1). The de- 
pendency on the exponent is basically the same as for A = 1 
and for sufficiently steep distributions the ejected mass frac- 
tion removes sufficient energy so that the BHs can coalesce. 
At 7 = 3 we can approximate Of/o g with 1 + In A/ In r\, a 
function which increases slowly with A. Even for the cluster 
extending as far as r c — rja^ only about twice as much mass 
is needed for the BHs to coalesce as compared to A = 1. 

To obtain the mass required for coalescence of the BHs 
for other exponents than 7 = 3, where exactly m c j is needed 
(see Fig. 1), we can proceed as follows. In Eq. (12) we assume 
the cluster mass to be distributed between a g and a h , i.e. 
A = 1 and C, — rj. Scaling the required mass to m c j with a 
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factor x = M(a h )/m c j we just need to replace 1/fc with \/k 
in Eq. (12) and solve ah/at = rj for x- This yields 

M(g h ) _ fc 2- 7 1-g 7 " 3 V-l , > 

m ej 3 - 7 1 - 7?t- 2 Inry 

Plotting this ratio as function of the exponent for k = 1 
again shows that for flat distributions (7 < 2) much more 
mass is required than for steeper cusps, at least for r\ > 
50 (Fig. 2). For 7 « 2.5 this is less than 2m e j for all ?j, 
confirming our conclusion that the binary is likely to merge 
in compact cusps which are steeper than 7 = 2. 

Steep density distributions with 7 > 2 have not been 
explored in numerical simulations or scattering experiments 
which used cusps at most as steep as 7 = 2. Considering the 
population of stars which is bound to the binary we obtain 
quite different results which show that the BHs can actually 
merge if the cusp is steep enough. 



BH, and is also in good agreement with fits to the mass 
deficit derived from observations (Milosavljevic et al. 2002). 
Because the mass bound to the binary is more than twice 
as much as the mass required for coalescence we conclude 
that the fraction of bound stars in the loss-cone contributes 
at least a significant fraction to the coalescence and might 
be enough to enable merging until a — a g on its own. We 
assumed that this mass is moving on circular orbits and ap- 
proximated the potential of the binary by that of a point 
mass. In Figs. 1 and 2 we showed that the ejection of mass, 
which is distributed in the potential of the binary according 
to a flat profile (7 < 2) as in previous papers, does not allow 
the BHs to merge. Either the binary stalls or a very large 
amount of mass is required. However, these figures also show 
that for steeper distributions the ejected mass is indeed suf- 
ficient to remove enough energy from the binary so that the 
BHs can coalesce. 



4 DISCUSSION AND CONCLUSIONS 

Although numerical simulations showed that the BHs in a 
hard binary coalesce after about 10 ~ 8 yr if the loss-cone is 
full (e.g., Quinlan 1996; Milosavljevic & Merritt 2001; Zier 
& Biermann 2001) most authors argued that the loss-cone 
becomes depleted and the binary stalls. This conclusion is 
based on the assumption of a spherically symmetric shallow 
density profile and has been confirmed recently by Berczik 
et al. (2005), using initial conditions which are in favour 
of a stalled binary. The choice of a flat central density dis- 
tribution is based on profiles derived from observations of 
elliptical galaxies, i.e. galaxies before or (more probably) af- 
ter a merger, when mass has been redistributed from the 
inner to the outer parts of the cluster, resulting in a flat- 
ter profile. However, during the merger the profile might be 
quite different and much steeper. Simulations by Toomre 
& Toomre (1972) support this conjecture and Barnes & 
Hernquist (1996) showed that both cores accumulate a big 
amount of mass before they merge and the binary becomes 
hard. Each of the BHs will carry a stellar cusp with a mass 
of about its own. 

Hence in this Letter we suggest that by the time the 
BHs become hard, mass and angular momentum have been 
redistributed and energy dissipated so that the central cusp 
is more massive and has a steeper density profile than that 
deduced from observations. Naturally this is a transient dis- 
tribution and unlikely to be observed because of its short 
lifetime. Because a mass of about 2 M12 will bind to the bi- 
nary, a large fraction of which will probably be contained 
in the loss-cone, we focused on this population, unlike nu- 
merical experiments which so far have focussed on scattered 
unbound stars. By computing the binding energy of a stel- 
lar distribution in the potential of a binary we could show 
that the ejection of this population can play a decisive role 
for a successful merger. From the energy extracted from the 
binary by these stars due to sling-shot ejection we derived a 
new expression for the mass which is required to be ejected 
so that the BHs coalesce (Eq. [8]). While previous work as- 
sumed this mass to be proportional to the mass of the bi- 
nary M12 due to the choice of the scaling factor, we find it 
to be proportional to the reduced mass. This seems to be 
more plausible, especially in the limit of a small secondary 



If we allow the stars to move on eccentric orbits us- 
ing their former radius as pericenter, the binding energy 
will decrease while the required ejected mass will increase. 
However, this will be more than compensated by the non- 
vanishing velocity of ejected stars at infinity, which we as- 
sumed to be 0. Comparing Eq. (6) with Yu (2002) we find 
k — 3.2, including elliptic orbits. Elliptic orbits allow the 
matter to be distributed less steeply than what is suggested 
by our results and therefore relax the conditions for a com- 
plete merger. We also assumed the BHs to move on circular 
orbits. If they are elliptic gravitational radiation dominates 
the shrinking earlier, shifting a g to larger radii and conse- 
quently reducing the mass required for ejection. Therefore 
our results in fact are quite conservative and we conclude 
that the BHs coalesce in most cases, contrary to previous 
results which were based on flat cores. Our analyses can also 
be applied to dark matter. If there is not sufficient baryonic 
matter for a merger and the BHs still merge, our approach 
allows some conclusions to be drawn about the central dis- 
tribution of dark matter. 

Simulations of merging BHs in a triaxial potential (Yu 
2002; Holley-Bockelmann et al. 2002; Holley-Bockelmann & 
Sigurdsson 2006; Berczik et al. 2006) strongly support the 
conclusion that the BHs merge completely in less than a 
Hubble time. Due to the triaxiality of the potential, the 
central regions remain populated. The flow of orbits through 
phase space into the loss-cone is large enough to keep it al- 
ways filled. An unfilled loss-cone is only found for shallow 
spherically symmetric cores (Yu 2002). This is an idealized 
distribution which, if at all, applies only to the cluster af- 
ter the merger has been completed, but not to an ongoing 
merger. Thus, additionally to the population bound to the 
binary that we were investigating here, simulations of a tri- 
axial potential show that the loss-cone remains filled during 
the merger, supplying the binary with an almost unlimited 
amount of stars that can be scattered off. Even in the case 
that the bound stars alone are not enough to allow the BHs 
to merge, together with the scattered stars from the loss- 
cone a stalled binary should be most unlikely. 

In a future paper (in preparation) we will discuss this 
approach and its context in more detail. 
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